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We study statistical properties of a one dimensional infinite system of coalescing particles. Each 
particle moves with constant velocity ±v towards its closest neighbor and merges with it upon 
collision. We propose a mean-field theory that confirms a t~^ concentration decay obtained in 
simulations and provides qualitative description for the densities of growing, constant, and shrinking 
inter-particle gaps. 

PACS numbers: 02.50-r, 01.75-fm, 89.90-|-n 



We introduce a simple deterministic model describing a 
coarsening dynamics of interacting domains. We arrived 
at this model in attempt to model the essential features of 
development of countries and civilizations. In our model, 
the civilizations are represented by domains that contin- 
uously cover one-dimensional "world" without gaps and 
overlaps. Neighboring civilizations are engaged in a per- 
manent warfare: A bigger civilization invades the lesser 
one, so that the border (interface) moves with velocity 
±t;/2 with V being the same for all pairs of neighbors. 
The model is equivalent to an infinite particle aggregat- 
ing system in which each particle (corresponding to the 
interface between domains) moves ballistically towards 
its closest neighbor. When a civilization shrinks to zero 
size and disappears, two particles, corresponding to its 
borders, collide. Such collisions between interfaces lead 
to coalescence, I + I ^ I where / symbolizes an inter- 
face, and the emerging interface starts moving towards 
its nearest neighbor with velocity of the same magnitude 
v/2 (Fig. 1). 
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Fig. 1. Schematic illustration of domain structure. Domains 
that grow, remain constant, and shrink at early time are de- 
noted by G, N, and S, respectively. 

The gradual character and general nature of the model 
makes questionable its ability to reproduce specific de- 
tails of human history that have been actually observed. 
However, it does provide some statistical characteristics 
like size distribution and predicts decrease of the total 
number of civilizations. An appealing simplicity of our 
model suggests that it might be relevant to description 
of other coarsening phenomena. 

In our previous work Q, we studied a model with 
domains growing freely in otherwise empty space and 



engaging in instantaneous warfare upon contact. One 
can view the model of Ref. |jl|] and the present model 
as two limiting cases of the general growth- and- war pro- 
cess: The former growth- controlled process is limited by 
the growth rate while the latter process is war- controlled. 
Another model which has strong relevance to our sys- 
tem is the "cut- in- two" model of Ref. ||], introduced 
in the context of breath figures coarsening. Similarly 
to our model, it assumes the complete coverage without 
overlaps, but unlike our model, coarsening events are in- 
stantaneous. Namely, evolution proceeds by consecutive 
elimination of the current shortest domain. 

In our model, which can be also called the ballistic 
coalescence model, there are three types of domains - 
growing, neutral (with constant length), and shrinking. 
A growing domain is bigger than both of its neighbors, 
so its length grows with velocity v. A neutral domain 
has one neighbor which is bigger and the other which is 
smaller, so its length remains constant though its location 
does change. Finally, a shrinking domain is smaller than 
both of its neighbors, and its length decreases in time 
with rate v. When such shrinking domain disappears, a 
change takes place in one of the two neighboring domains 
- a growing domain may become neutral, or a neutral do- 
main may turn into a shrinking one. The number density 
of growing domains, G(t), is always equal the number 
density of shrinking domains, S{t), since the space re- 
mains continuously covered without gaps and overlaps 
throughout the process. A growing domain of initial 
length £q has the length vt -\- £o at time i, implying that 
G{t) ^ and S{t) ^ t~^. One might expect that the 
number density of neutral domains, N(t), decays slower 
than t~^, implying that the total domain density also 
exhibits anomalously slow decay. This indeed happens 
e.g. in apparently similar model of ballistic annihilation, 
/ -I- / 0, where t^^/^ decay of the domain density is 
observed . In both ballistic annihilation and ballistic 
coalescence models, growing domains cannot be neigh- 
bors, while an arbitrary number of consecutive neutral 
domains may coexist. However, in the ballistic annihila- 
tion model such a train of neutral domains has a good 
chance to live long while in our model this sequence of 
neutral domains can be eliminated by a single growing 
domain. Because of this relative vulnerability of neutral 
domains in the ballistic coalescence model, their num- 



1 



ber density scales similarly to the other two densities, 
G{t) and S{t). Moreover, we will show below that N{t) 
becomes smaller than G{t) and S{t) as the process de- 
velops. Note that initially all the densities are the same, 
G(0) — N{0) — S{Q). Indeed, the initial sizes of domains 
are uncorrelated, so if we take three consecutive intervals 
the probabilities that the middle interval is the smallest, 
medium, or biggest are all the same. 

Thus in the long time limit the densities behave ac- 
cording to 
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with constants g — s and n to be determined. Con- 
sider an arbitrary shrinking domain of vanishing length. 
It may be surrounded by two neutral, two growing, or 
growing and neutral domains. Denote by A, /i, and 
the respective probabilities {X + fi + ly — 1). The rate 
at which domains disappear, vS{0,t), may be written as 
ag/t^ with some constant a. Here S{x,t) is the den- 
sity of shrinking domains of length x. The normalization 
/p dxS{x,t) = S{t) - t-^ suggests S{Q,t) - t'"^ and ex- 
plains the rate given above. Now it is straightforward to 
write down the rate equations 
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Summing Eqs. (g) gives a = {g + n + s)/g = 2 + n/g. The 
above equations are exact. To proceed further we need to 
know A, /i, ly. We assume that the types of domains sur- 
rounding vanishing shrinking ones are uncorrelated, so 
the probability of picking up growing or neutral domain 
is proportional to its respective concentration. Mathe- 
matically, it means that 
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Eqs. (g) imply = 1 and a{l 
ing a = 2 + n/g and fj, = g'^{g - 
equations allows us to find n/g, 



- 2/i) = n/g. Substitut- 
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n 
9 



V5 - 1 



P, 



(4) 



where p is known as the "golden ratio". We still need 
to determine one constant, say g. This can be accom- 
plished by using the fact that the system is continuously 
covered with domains. We need to determine the length 
distributions G{x,t), N{x,t), and S{x,t). In the limit 
of large time, we may ignore differences in sizes of grow- 
ing domains. These differences are determined by initial 
size distribution and do not change with time. Thus the 
length distribution of growing domains is simply 



G{x,t) ^^5{x-t). 
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Two other length distributions are expected to scale, so 
we write 



iVKO = |^(i), S(^,«) = |B(f) 
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In the mean- field approximation (Eq. (^)), the govern- 
ing equations for N(x, t) and S{x, t) read 
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The spatial derivative term in (^) accounts for continu- 
ous shrinking. The velocity of shrinking v is set equal 
to one without loss of generality. We also have dropped 
two terms proportional to 5-functions - the gain term in 
Eq. relevant for conversion of growing domains into 
neutral, and the loss term in Eq. (@) accounting for re- 
moval of zero-size shrinking domains. These terms pro- 
vide appropriate boundary conditions. 

Using the scaling form (^) and relations (|^) we simplify 
(0) and i) to 



uA' ^2A dvA{v) 



and 



2B + {1 + u)B' = -2A ~2A dvA{v) 
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In these equations u = x/t, A' ^ dA/du, B' = dB/du. 
According to the definitions (|l|) and @ , the scaling func- 
tions A{u) and B{u) obey normalization conditions 



duA{u) — p, 



[ duB[u) = l. 
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Additionally, creation of the longest neutral domains 
(x = t) and removal of the zero-size shrinking domains 
give N{t,t) = apg/t^ and S'(0,i) = ag/t^. This provides 
the boimdary conditions: 



A{l) = l, B{0) = a. 



Solving (0) gives 



A{u) = 
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Note that A{0) = 0, in agreement with intuitive expec- 
tation that there are no neutral domains of length zero. 
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To solve for the length distribution of shrinking do- 
mains we first multiply (H) on 1 + u so that the left-hand 



side becomes a complete derivative, [(1 



The 

right-hand side contains already known functions. Inte- 
grating the resulting equation with the boundary condi- 
tion B{0) — a = 2 + p we obtain 
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One can check that B{1) = 0, the result that agrees 
with intuitive expectation that there are no shrinking 
domains of maximal possible length x = t. We verified 
that the normalization condition of Eq. ( pi] ) is satisfied. 
The scaled length distributions of neutral and shrinking 
domains are plotted on Figs. 2 and 3. 

Now we can compute the constant g. The requirement 
that the space is completely covered, 



dxx[G{x, t) + N{x, t) + S{x, t)]^l, 



(15) 



gives 



duu[A{u) + B{u)] 



0.625837. (16) 



We performed molecular dynamic simulations for sys- 
tems of 2 • 10^ particles with periodic boundary condi- 
tions. Results were averaged over 500 initial configura- 
tions. We found that the ratio n/g of neutral to grow- 
ing particle densities decreases slowly from 1 to approxi- 
mately 0.66. This value is 7% larger than the theoretical 
prediction n/g = p ^ 0.618033989; the statistical error 
for the measured value is less than 3%. 
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Fig. 2. Plot of scaled domain size distribution for neutral 
domains, A{x/t). Time is incremented by 2.25: t = 2192(0)j 
t = 4932(o), t = 11097(+), t = 24968(*). Mean-field predic- 
tion (Eq. 13) is shown by dashed line. 



We also determined numerically the probabilities 
A, /i, J' and found more pronounced difference with mean- 
field predictions of Eq. (|): A w 0.1, p, « 0.5, v « 0.4 

vs. Amft = 0.145898034, ^MFT = 0.381966011, j^mft = 
0.472135955. We found that scaling (||) for the neutral 
and shrinking domain length distribution functions holds 
indeed, although the shape of the experimental curves is 
different from mean-field predictions (Figs. 2 and 3). The 
discrepancy is particularly drastic for the size distribu- 
tion of neutral domains indicating that correlations, not 
accounted by the mean-field theory, are especially impor- 
tant for them. 
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Fig. 3. Plot of scaled domain size distribution for shrinking 
domains, B{x/t). Time is incremented by 2.25: t — 2192(0)) 
t = 4932(o), t = 11097(+), t = 24968(*j. Mean-field predic- 
tion (Eq. 14) is shown by dashed line. 
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Fig. 4. Schematic illustration of domain structure for a 
system with interface annihilation 

A natural generalization of this system is a model of 
an "alternating bipolar world" . Here neighboring civi- 
lizations are hostile to each other so they are engaged in 
warfare, while collisions of friendly civilizations lead to 
merge. In other words, interfaces annihilate upon con- 
tact rather than coalesce as illustrated on Fig. 4. When 
a shrinking domain disappears, its left-hand side and 
right-hand side neighbors merge into a single domain. 
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Since a domain length can instantaneously increase upon 
each shrinking event, such events often lead to changes 
in as many as four domains surrounding the shrinking 
one (central group of domains on Fig. 4). The model is 
thus different from the ballistic annihilation model 
Furthermore, correlations become more important than 
in our original ballistic coalescence model. However, this 
interface annihilation model exhibits the same scaling as 
the original interface coalescence model, i.e., the length 
of a typical growing, neutral, or shrinking domain grows 
linearly in time, and domain concentration decays as 1/t. 
Scaled plots of length distribution for growing, neutral, 
and shrinking domains are presented of Fig. 5. Note that 
unlike the original model, the length distribution of grow- 
ing domain remains unsingular and unbounded. 
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Coarsening rules in the above models are similar 
to those in one-dimensional Potts models with zero- 
temperature Glauber dynamics. Indeed, in the infinite- 
state Potts model interfaces aggregate upon collisions 
while in the 2-state Potts model, i.e. the Ising model, 
interfaces annihilate upon collisions. The difference lies 
in interface dynamics - in the Potts models domain walls 
usually undergo random walk rather than ballistic mo- 
tion. 

In summary, we have analyzed a model of ballistic ag- 
gregation of particles where each particle is approaching 
its nearest neighbor with constant and universal veloc- 
ity. It can be considered as an idealized model of a world 
in which all neighboring countries are engaged in perma- 
nent warfare with larger countries advancing and smaller 
receding. We found a universal linear in time scaling 
for all types of domains and developed a mean-field like 
approach to the domain size distribution function. 
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Fig. 5. Plots of scaled domain size distribution for grow- 
ing domains (o), neutral domains (0)i S'Hd shrinking domains 
for interface annihilation model. Simulation included av- 
eraging over 500 configurations of 2 million initial domains 
each. 
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